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Abstract 

We perform two independent calculations of the two-loop partition function for the 
large N 't Hooft limit of the plane-wave matrix model, conjectured to be dual to the 
decoupled little string theory of a single spherical type IIA NS5-brane. The first is via a 
direct two-loop path-integral calculation in the matrix model, while the second employs the 
one-loop dilatation operator of four-dimensional M = 4 Yang-Mills theory truncated to the 
SU(2\4) subsector. We find precise agreement between the results of the two calculations. 
Various polynomials appearing in the result have rather special properties, possibly related 
to the large symmetry algebra of the theory or to integrability. 



1. Introduction 

A recent and fascinating addition to the cast of maximally supersymmetric Yang-Mills 
theories has been the plane- wave matrix model JIJ . This theory, a massive deformation of 
the BFSS Matrix Theory preserving all 32 supercharges, has been conjectured to describe 
(in a particular large N limit) M-theory on the maximally supersymmetric plane- wave so- 
lution of eleven-dimensional supergravity. The theory turns out to be much more tractable 
than the usual BFSS matrix model p|, allowing perturbative computations at fixed N in 
the limit of large mass m and harboring a powerful symmetry algebra [|3],|],[5| that allows 
extrapolation of some perturbative results into the strongly-coupled regime. Among this 
reliable information at strong coupling is direct evidence that certain vacuum states of the 
model describe spherical BPS transverse M5-branes of M-theory ||. 

Recently [|7|,§, there has emerged an interesting connection between the plane- wave 
matrix model and four-dimensional Af = 4 Super- Yang-Mills (SYM) theory on S s . At 
the classical level, the plane-wave matrix model emerges directly from the Af = 4 theory 
through a consistent truncation |7| that keeps only the modes invariant under a certain 
SU(2) subgroup of the SU(2, 2|4) superconformal algebra (the "S77(2|4) subsector"). In 
the strict large N 't Hooft limit, the relationship between these two theories extends further. 
The set of states of the Af = 4 theory built from these modes forms a subsector which 
is closed under renormalization at one loop. Surprisingly, all one-loop corrections to the 
energies of these states in the Af = 4 theory match precisely with the one-loop correction 
to the energies computed in the plane- wave matrix model after the correct identification 
of couplings between the two theories |7|,[| . The common one-loop Hamiltonian governing 
these energy shifts corresponds to an integrable SU{2\4) spin chain, so at least at the 
one-loop level, the recently discovered integrability properties of the planar Af = 4 SYM 
theory extend to the plane-wave matrix model. In fact, recent explicit calculations §| in 
the plane-wave matrix model (for a certain closed subsector of scalar modes) suggest that 
both the integrability properties and the equivalence with the appropriate subsector (the 
U SU(2) subsector") of the Af = 4 theory persist even to three loops! 

Given these results, it is natural to wonder whether the full plane-wave matrix model 
is integrable in the 't Hooft large N limit. It is important to note that this limit is quite 
different from the Matrix Theory limit conjectured to define M-theory on the plane wave 
(which clearly has no chance of being integrable). In fact, it was argued in || that the 't 
Hooft large N limit, defined about the trivial vacuum state of the theory, is a decoupling 
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limit which only keeps the excitations of a single spherical fivebrane. More precisely, since 
this limit does not decompactify the M-theory circle, the fivebrane should be interpreted 
as a spherical type IIA NS5-brane with the 't Hooft parameter related to the sphere radius 
in units of a'. The fact that this limit describes the decoupled physics of only a single 
brane (according to the arguments of 0) provides additional hope that it may indeed be 
integrable. 

If the 't Hooft large iV limit of the plane-wave matrix model does turn out to be 
integrable, one might aspire to calculate the exact spectrum for all values of the coupling, 
or equivalently, to find an analytic expression for the exact partition function as a function 
of coupling. Motivated by the hope that such an expression exists, we proceed in this note 
to calculate directly the leading terms in its weak coupling expansion. Thus, we compute 
the two-loop partition function for the strict large iV limit of the plane- wave matrix model 
about its trivial vacuum,0 extending the previously calculated fK||llJ zero-coupling result 
given by (2.6) below. We find that the correction to the partition function takes the form 

41 U=1 l-z((-l)"+ V) q(y)f 

where A = ^m-^ ^ s the 't Hooft parameter of the Af = 4 gauge theory (we relate it to 
the matrix model parameters below), y = e~ /3rn ^ 12 , and z, g, p and q are polynomials in y 
given below in (2.5), (2.20) and (2.24). These have some rather special properties that we 
comment on in section 4. 

Our calculation is carried out by two independent methods. The first method, in 
section 2, is a direct two- loop path-integral calculation using the Euclidean matrix model 
action with Euclidean time compactified on a circle of radius j3 = 1/T. Our second method, 
described in section 3, amounts to an explicit sum over states of the Boltzmann factor, 
taking into account the leading order energies together with their one-loop corrections. 
For this approach, we use the one-loop equivalence to the SU (2|4) subsector of the N = 4 
SYM theory, and apply the general analysis of to express the subleading terms in 
the partition function in terms of the one-loop dilatation operator of four-dimensional 



1 At N = oo with fixed finite 't Hooft coupling, the free energy diverges at a finite temperature, 
and our results for the partition function are valid below this temperature. At large but finite N, 
the story is more complicated, since the model has of order vacua which should all contribute 
since their ground state energies are all zero. Thus, our result should only be interpreted in the 
context of the strict large N 't Hooft limit for which the vacua decouple. 
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M = 4 Yang-Mills theory truncated to the SU (2|4) subsector. While the details of the two 
calculations look rather different, both calculations precisely give ( |1 . 1| ) . 

Even in the limit of zero coupling (studied previously in |jIl,lQ|j), this partition function 
displays interesting Hagedorn behavior, with a limiting temperature in the strict large 
N limit at which the free energy diverges logarithmicallyi This Hagedorn behavior is 
presumably associated with the Little Strings of the decoupled IIA NS5-brane defined 
by this limit. From our two-loop results, we can determine the change in the Hagedorn 
temperature as the coupling is turned on (i.e. as the sphere on which the Little Strings 
live grows from zero size), and we find that it increases with the coupling for small A. This 
is consistent with the suggestion [|J that the strong coupling limit should be equivalent to 
the free conformal theory associated with a single spherical M5-brane, for which we expect 
no Hagedorn behavior. 



2. Plane-Wave Matrix Model Path Integral 

In this section, we compute the two-loop partition function for the 't Hooft limit 
of the plane-wave matrix model directly via a path-integral calculation. We follow all 
matrix model conventions of ||, in which the plane-wave matrix model action in Euclidean 
signature is given byi 

C = Tr QpoXT + ^(X? + \{D G X«f + ^(X«f + i^ Ia D ^ Ia + j^ Ic }) 
+ R-*Tr (^e ijk X i X^X k + ^ Ia a i J[X i ^ I/3 ] 

~e a p^ aI gh[X a ,^ J ] + l -e a ^ aI g^ IJ [X a ^ aJ ]^ 
- R 3 Tr (^[X\X J ] 2 + ^[X a ,X b ] 2 + ^[X\X a } 2 ^j . 

(2.1) 



For large but finite N, this divergence signals a phase transition to a deconfined phase with 
free energy of order TV 2 [13,11,10|]. 

3 Scalar indices k and a,b,c are associated with the vector representations of SO(3) and 
50(6) respectively, while fermion indices of a, (3 and /, J are in the spinor representations of these 
groups. The g's are Clebsch-Gordon coefficients relating the vector of S0(6) to the antisymmetric 
product of two spinors (fundamentals of SU(4)). We set lp = 1, but we can restore lp in any 
formulae using the fact that R and 1/m have dimensions of length. 
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Choosing the gauge d t A Q = 0, and introducing the corresponding Fadeev-Popov de- 
terminant A, the thermal partition function is 

Z = J[dX l ][dX a }[dip Ia }[dA }Ae~ So dtc ^ (2.2) 

where the time direction has been compactified with radius (3 = 1/T and bosons/fermions 
are taken to have periodic/antiperiodic boundary conditions respectively. As explained in 



Tl],|T0| , all modes are effectively very massive at weak coupling except the zero mode of the 
gauge field on the thermal circle. It is then convenient as an intermediate step to integrate 
out all other modes to produce an effective action for this zero-mode, which we denote 
by a. As argued in [II]], the resulting effective action for this mode may depend only 
on the unitary matrix U = e ll3a (the Wilson line of the gauge field around the thermal 
circle), and the effect of the determinant A in the path integral is precisely to convert 
the measure [dAo] into the Haar measure [dU] for unitary matrices. Thus, the partition 
function reduces to an ordinary unitary one-matrix model. 

2.1. One-Loop Result 

The evaluation of the partition function to one 1-loop has been carried out in [|TT|,|TU 



The result is 



where 



Z x - 



p = J [dU}e- s °->°°*W (2.3) 
^i-ioop(^) = E -^((-l) n+ V)Tr(^)Tr(^). ( 2 .4) 



n 

n=l 



Here, we define y = e P™-/ 12 and 

z(y) = 6y 2 + 8y 3 + 3y 4 (2.5) 

is the single mode (letter) partition function. At strictly infinite N, the free energy has a 
Hagedorn divergence at Th = m/(121n3). Below this temperature, the model is governed 
by a stable saddle point for which the eigenvalues of U are spread uniformly around the 
unit circle, so that Tr(£7 n ) = 0. Performing the Gaussian integral about this configuration 
gives the first non-zero contribution to the free energy, and we obtain 

OO j 

z 1 -.^=«- w --=n 1 - z(( -i ) ^ V) - < 2 - 6 > 



2.2. Two-Loop Calculation 

At two loops, the partition function is given by 



^2-ioops = J [dU]e- s °-ioo P (u)-s?l loop (u) > (2 _ 7) 



where 

^2-loop 



g °2-loop = /g °mt \ 



(2.* 

exp(-(^ 4 ) + i^3^3)), 



where £3 and £4 are the cubic and quartic terms in the action ([2J]). Here, the expectation 
values are evaluated in the free theory with fixed a. 

The required propagators follow from the quadratic action 



S 2 = Jdt Tr {^X\-Dl + ^-)X* + \X«(-Dl + + ^(D + |)^«) . 

(2.9) 

For the boson propagators, we find 

2.10 

(X; q (t)X b rs (t')) = S ab A f (t - t>, a)^. 
The propagator A is defined to be a periodic function of t given in the domain [0, f3) by 

e iat / e~ Mt e Mt \ 

A M (t, a) = — I — — j-^Mp - — — w^MpJi ( 2 - n ) 

where a is short for (a ® 1) — (1 ® a) and matrix indices have been suppressed. The fermion 
propagator is 

((V'/ a (t)) P9 (^ tJ/3 )(t , ))=^fA|(t-t',ar^, (2.12) 
where A F is defined to be an antiperiodic function of t given in the domain [0, (5) by 

-Mt 

A M(^) = e» t T --^. (2.13) 
There are six correlators contributing to £2- loop- These are 

= J dtTr([X\X^)), (2.14) 

J dtTr([X\X a } 2 )), (2.15) 



S2 = 7T 



/ ^ 3 



;2.i6) 



J dtTr([X a ,X b } 2 )} 
s 4 = (^^- J dt Tr(e ljk X i X j X k ) J dt' Tr(e lmn X l X m X n )) 

55 = ( ~1T / dtTr^a'JiX',^}) J dt' T^ tJ ' Q </ [X\ ^,\)) 

56 = ( l~/ dt ^M ]Ia gh[X a A UP ]) J dt' Tr(e a '^^ Ila ^ a ' I,J '[X a ' 'rfj.y])). 

These correlators contribute to the partition function in two different ways. First, 
planar diagrams contribute terms to the double-trace effective action for U which modify 
the Gaussian integral and result in order A corrections to the denominators in (|2.6|) . Since 
the Gaussian integral is actually the subleading contribution to the large iV free energy 
(the leading 0{N 2 ) contribution vanishes for this saddle point) there are also contributions 
at the same order arising from nonplanar two-loop diagrams.^ These are independent of 
U and give a temperature-dependent prefactor to the infinite product in ( |2.6| ). 

2.3. Planar Contribution 

From the six correlators above, we first write the planar contributions, giving in 
the first line the complete expression for the planar part of the correlator in terms of 
propagators and in the second line, the terms contributing to the double trace action. 
There are in addition three-trace terms, but these do not contribute to the partition 
function at infinite N. We find: 



(ai)pi = 3(3R 3 A^ (0, a ab ) Am (0, a ac 



?3 

T 

oo 



27 Ft^ AT 

V —p — {y Sn + 2y^)Tr(U n )Tr(U^) 

< * ) rn z 



, 2 m 

n=l 

(a 2 ) P i = 18/3i? 3 A f (0,a ab )A f (0,a a 

oo 



V 162/3^-^ (y 6n + y 4n + y 2n )Tr(U n )Tr(U^ n ) ' ~ ' ' ' 

' J qnrt — 



m 

n=l 



s 3 ) p i = 15/3i? 3 A t (0,a a6 )A^(0,a ac ) 



oo ^3 



Ft IV 

V 270p—^ r (y 4n + 2y 2n )Tr{U n )Tr(U^ 

n=l 



We thank Ofer Aharony for emphasising that non-planar diagrams must play a role here. 



(s 4 ) P i = -m 2 [3R 3 J dtAm_(t, a a6 ) A f (£, a bc ) A f (t, a c 

- E t^^tV 71 - % 4n )^(^ n )Tr(^) 



2 m 

n=l 

(s 5 ) p i = -24/3i? 3 y dt A - (t, a a6 ) A| (t, a 6c ) A| (/3 - t, a ac ) 

00 D3/Y ( 2 - 18 ) 

-> E 216/3^- (-l)^ 1 ^ 6 " + 2y 5n + 2y 3n - 3y 2n )Tr(t/ n )Tr(t/ tn ) 

n=l 

(s 6 ) P i = -48/3i? 3 J dtA f (t, a ab ) A| (£, a 6c ) A| (t, a ca ) 

-> V 216/3^4- (-l) n+1 (y 7n - y 3n )Tr(t/ n )Tr(£/t n ). 

n=l 

In the expressions above, each of the propagators contribute factors of a to two of the three 
index loops, which we label by a, 6, and c. The notation a ab indicates that for the tensor 
products (a <g) 1) — (1 <g) a) appearing in the propagator, the first and second elements of 
the tensor product appear in the traces associated with index loops a and b respectively. 

Combining all terms, we find that the two-loop contribution to the double-trace effec- 
tive action for U is 

00 

S?: ioop (U) = -\J2 My)g((-iy +1 y n WU n )Tr(U^), (2.19) 

where 

g(y)=y 2 (l + y) 4 (l + y 2 ) (2.20) 
and we have defined a 't Hooft coupling A = 648R s N/m 3 . 

2.4- Nonplanar Contribution 

We now evaluate the nonplanar contributions from the six correlators above. In 
this case, since there is only a single index loop and since each term in the propagators 
contributes an equal number of Us and U^s, we will always end up with just the identity 
matrix inside the single trace. Thus, the same result will be obtained by setting a = 
(U = 1) in all propagators from the start. The nonplanar contributions are thus 



(si)np = -3^ 3 A f (0,0) A™ (0,0) 



27 R 3 N (1 + y 4 ) 2 (2-21) 
4 1 m 2 (1-y 4 ) 2 ' 
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(s 2 ) np = -18(3R 3 A f (0, 0) Am (0, 0) 
ff 3 iV(l + y 4 )(l + y 2 ) 

61P W 2 (l-y4)(l- y 2) 

(s 3 )n P = -15/3# 3 A ? (0,0)A f (0,0) 
,i? 3 iV(l + y 2 ) 2 



- 135/3- 



m 2 (1 — y 2 ) 2 



/np — "« H*-* I "<>"-* - 

>3 /\r n _i_ m„4 



(s 4 ) np = m 2 /3i? 3 / M f (t,0)A f (t,0)A f (t,0) 



- 27 i? 3 iV(l + 10y 4 + y 8 ) (2.22) 
4 tm 2 (l-y 4 ) 2 

(s 5 ) np = 24(3R 3 y eftA f (£, 0) A| (£, 0) A| (/3 - £, 0) 

~ bP m 2 (l + y 3 ) 2 
(se)np = 48/3.R 3 ^ dtA™ (t, 0) A| (£, 0) A| (£, 0) 

m 2 (1 + y 3 ) 2 (l - y 2 ) 
Combining all contributions, we find 

S np = X\n(y)^ (2.23) 



where A is defined as above and 



p(y) = y 2 (l + y + y 2 )(l + y + 6y 2 + y 3 + 6y 4 + y 5 + y 6 ), 
q(y) = (1 - y) 2 (l + y) 2 (l - y + y 2 ) 2 (l + y 2 ) 2 . 



2.5. Summary: Two-Loop Result 



(2.24) 



Using the results above, it is now straightforward to complete the calculation of the 
two-loop partition function by performing the Gaussian integral around the Tr(£/" n ) = 
saddle point. Combining the one- and two-loop effective actions for U, we find that the 
terms quadratic in traces are 

Seff>) = E - f 1 " ^((-l) n+ V) - Aln(^)<7((-ir+V) + 0(X 2 )) Tr(U n )Tr(U^). 

n—1 

(2.25) 



From this, we can read off the appropriate modification of the denominators in (|2.6|) , 
so combining the results of the Gaussian integral with the prefactor coming from the 
nonplanar diagrams, we obtain our final result 



~- oo 
-X\n(y)p(y)/q(y) 



z((-l) n+1 y n ) - \ ]n(y n )g((-l) n + 1 y n ) 



+ 0(A 2 ), (2.26) 



where g, p, and q were defined in (|2.20| ) and fl2.24| ). For what follows it will be convenient 
to write the first O(X) correction as 

3. Dilatation Operator in the SU (2\4) Subsector 



In this section we obtain the partition function (|2.27|) from a one-loop calculation in 
M = 4 SYM gauge theory on M x S 3 by making use of the one-loop isomorphism between 
the plane- wave matrix model and the SU (2|4) subsector of the gauge theory. This subsector 
consists of those operators built out of the six scalar fields <fi, the eight positive chirality 
spinors A, and the three self-dual components of the field strength tensor F. Henceforth 
we will use the matrix model notation in referring to these modes respectively as X a , ipi a: 
and X 1 . Covariant derivatives of these fields do not appear in this subsector. 

The free Hamiltonian Hq of the matrix model is identified with the tree-level dilatation 
operator D in this subsector of gauge theory according to 

Til 

H = ^D . (3.1) 
o 

Quantum corrections to the matrix model effective Hamiltonian can be computed for large 
m via ordinary degenerate quantum mechanical perturbation theory in the parameter 
1/m 3 . The first order correction has been shown to agree with the one- loop correction 
to the gauge theory dilatation operator, 

Ho + -^Vf = ^(d + ^D 2 ) , (3.2) 

where A = g\^\N is the gauge theory 't Hooft parameter. Although the SU (2|4) subsector 
of the gauge theory is not closed under renormalization beyond one loop, we can take 
H = (m/6)D in calculating the partition function 

Z = Trfe-^] = Tr[y 2D ] + 0(X 2 ) (3.3) 
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since at the moment we are not interested in the higher order terms. 

To calculate the leading term Tr[y 2i:)o ] in the partition function one simply has to 
enumerate the operators appearing in the SU{2\A) subsector weighted by their bare di- 
mension. A straightforward application of Polya theory yields the result (|2.6|). The first 



correction to Tr[y 2 ' Do ] was studied in [jT2], where a combinatorial analysis of the anomalous 



dimensions of gauge theory operators revealed that the result has the general structure 



2tt 2 ^ 



n=l 



n(D 2 ((-l) n+1 y n )) >r^, -v m+1 m , n n+l,.nN\ 



m=l 



(3.4) 



This formula is valid for temperatures below the Hagedorn temperature at which the 
N = oo free energy diverges, as evidenced by the appearence of a pole at y = 1/3 in 
the first term of ( |3.4]) . Intriguingly, it is already apparent that this expression has strong 
similarities with ( p.27|) . 

The quantities (D 2 ) and (PD 2 ) are defined as follows. The one-loop dilatation oper- 
ator D 2 only acts on two neighboring fields in any single-trace operator, 

L 

D 2 Tr[Ai ...A L ] = J2 " ' ' D 2 (AiA i+1 ) ■ ■ ■ A L }. (3.5) 

i=i 

In the SU(2\A) subsector, each letter A\ corresponds to one of the 17 fields {X a , ipi a , X % }, 
so we can think of D 2 as a 289 x 289 matrix acting on a pair of letters IA1A2). The 
ingredients appearing in (|3.4|) are just traces of this matrix, 



(D 2 (y)) = Tr[y 2D °w +2D °vW 2 }, (PD 2 (u,w)) = Tr[u 2D °^w 2D °M PD 2 ], (3.6) 

where P\AiA 2 ) = (—l) FlF2 \A 2 Ai) is the (graded) permutation operator and Z? («) gives 
the bare dimension of the zth letter. 

The one-loop dilatation operator D 2 in the SU (2\4) subsector takes the form 

D 2 \X a X b ) = -\X a X b ) - -\X b X a ) + -5 ab 5 cd \X c X d ), 
m 1 2 ' 1 2 1 1 4 1 1 

D 2 \^ Ia X l ) = \i> Ia X l ) + l -e^ k ^j(\tp I(3 X k ) + \X k ip If3 )), 

4 V 7 (3.7) 

D 2 \X^ Ia ) = |XVla) + y jk CT j a l3 (WipX k ) + |*Vl/3>), 

D 2 \X a X') = \\X a X l ) - \\X*X a ) + -g^a^faj^jp) - |^i«>), 
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(3.J 



D 2 \X l X a ) = \\X'X a ) - \\X a X*) + l -a^g^ IJ {\^ Ict ^j P ) - \^ Ia ) 
3 111 

D 2 \^Ia^Jp) = ^Ialpjp) + jllpjplpla) + ^ I f3^ J a) ~ j\^J a ^Ip) 

- Y 2 9 a iJ<?(\x a x l ) + \x l x a )), 

D 2 \X a ^ Ia ) = ^\X a ijj Ia ) - ^i a X a ) + ^g ab I J (\^ Ja X b ) + \X b ^ Ja ) 

D 2 \X l X j ) = -\X l X j ) + -\X J X l ) - -5 ij 5 kl \X k X l ). 
i ' 4 1 4 2 

We obtained these expressions by restricting the general result written in to the 577 (2|4) 

subsector. The first line, familiar as the Hamiltonian of the SO (6) spin chain, has been 

reproduced by a direct calculation in plane-wave matrix perturbation theory |7j . It would 

be interesting to extend their analysis to the full SU (2\4) subsector. The last line is the 

standard spin-1 SU(2) spin chain Hamiltonian. The required traces ( |3.6| ) can be easily 

read off from these formulas, and we find 

(D 2 (y)) = ly\l + y) 3 (11 + 72/), 

2 9 (3-9) 
(PD 2 (u, w)) = --u 2 w 2 (l + u)(l + w)(3 + u + w - uw). 

It remains to plug (|3.9| ) into ( |3.4[ ). The first trace can be broken into the two terms 

(D 2 (y)) = \g{y) - ?(1 - z(y))y 2 (l + 4y + 2y 2 ), (3.10) 

where g(y) is the same function we defined in ( |2.20| ). The factor of 1 — z in the second 
term cancels the denominator in (|3.4|) , allowing this term to be summed explicitly, giving 



v n(£> 2 ((-l)"+V)) _ 3 ^ ng((-l)"+y ) 

^ 1 - z((-l)"+V) 2 4^ 1 - *((-l) n+ V) ' 



n l ■ n l 

with 



rr x 3y 2 (l + y + r/ 2 )(l + y-3y 2 + 4r/ 3 -3r/ 4 + y 5 + y 6 ) 
U2/J 2(l-y) 2 (l + y) 2 (l+y 2 ) 2 (l-y + y 2 ) 2 ■ ^ > 



The double sum of the second trace in ( p^4T ) can be evaluated explicitly, giving the contri- 
bution 



9y 4 (l + y + y 2 )(3-y + 3y 

2(l-y)2(l + y)2(l +y 2)2(l_j / + y 5 
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Remarkably, Ti(y) + T 2 (y) = —p(y)/q(y) (defined in ( |2.24j )), so combining all terms yields 
the final expression 



SlnZ 



3Aln(y) J ^ ng((-l) n+1 y n ) p(y) 



, l - z((-i)™+V) g(y) I ' °' ' 



.n=l 



in precise agreement with (|2.27 ). The overall coefficient agrees after we make use of 
A = (47r 2 /3)A, which follows from the familiar relation 

/^3 = 32t^ 
\3R/ 9ym 

between the matrix model and Yang-Mills parameters. 



4. Discussion 

The result we have derived contains a significant amount of physical information about 
the model. Specifically, the coefficient of y n ln(y) in ( 2.26Q gives 12 jm times the sum of 



one-loop corrections to the energies of all states with energy nm/12 at zero coupling. These 
energy shifts result in a shift in the Hagedorn temperature at which the partition function 



diverges, and using (|2.27|) we find that the corrected Hagedorn temperature 

„ m ( 10 



H 



12 In 3 V 8l7r- 



i + ^A + OM , (4.1) 



increases as we move to stronger coupling. 

A particularly mysterious feature of our analysis is that both methods of calculation 
split naturally into two parts, yet the pieces on the two sides are not in direct corre- 
spondence. The two quantities g and p/q have quite distinct origins (from planar versus 
non-planar diagrams) in the calculation of section 2, while (-D2) and {PD2) have similarly 
distinct interpretations in section 3. However, the crucial equation ( |3.10| ) shows that there 
is not a direct identification between g and (D2). It would be very interesting to find a 
direct interpretation for g, or equivalently, to determine whether there is a sense in which 
the decomposition ( |3.10| ) is natural, from the spin chain viewpoint. 



Finally, it would be interesting to understand whether any features of the results we 
have derived are related to the special properties of the model, such as its large superalgebra 
or integrability properties. The overall form of the expression (|1.1|) is rather generic for 
large N gauged matrix models in 0+1 dimensions, so any special features should show up 



12 



in the polynomials g,p, and q themselves. In fact, these polynomials are quite non-generic: 
they possess a significant degree of factorization, a symmetry under reversing the order of 
exponents, and, in the case of g and q, all nonzero roots lie on the unit circle. Whether 



any of these features, some of which appear also in the full J\f = 4 SYM result fl^l , relate 
to integrability or supersymmetry is left as a question for future work. 
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